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We investigate polaron formation in a many-electron system in the presence of a local repulsion 
sufficiently strong to prevent local-bipolaron formation. Specifically, we consider a Hubbard-Holstein 
model of interacting electrons coupled to dispersionless phonons of frequency ujo- Numerically solving 
the model in a small one-dimensional cluster, we find that in the nearly adiabatic case uq < t, the 
necessary and sufficient condition for the polaronic regime to occur is that the energy gain in the 
atomic (i.e., extremely localized) regime Spoi overcomes the energy of the purely electronic system 
Sei- In the antiadiabatic case, ujo > t, polaron formation is instead driven by the condition of a large 
ionic displacement g/ujo > 1 (g being the electron-phonon coupfing). Dynamical properties of the 
model in the weak and moderately strong coupling regimes are also analyzed. 
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I. INTRODUCTION 

There is increasing evidence that polarons are formed 
in the lightly doped insulating phase of the high tem- 
perature superconducting cuprates ||l|,^ and in the high- 
temperature paramagnetic phase in doped manganites 
P|. On the theoretical side, simple polaronic models like 
the Holstein and the Su-Schrieffer-Heeger models have 
attracted intensive studies in the last few years. How- 
ever, whereas the cases of one or two polarons have been 
carefully investigated both numerically P-pO|| and ana- 
lytically ||ll|,|l^, the case of many polaronic carriers is 
still incompletely understood. The relevance of electronic 
correlations in the materials mentioned above provides a 
strong motivation to study the combined effect of the 
electron-phonon coupling which leads to polaronic fea- 
tures and the strong electron-electron interaction. The 
main goal of the present work is to numerically investi- 
gate strongly correlated many-particle systems interact- 
ing with phonons in order both to identify a simple crite- 
rion for the formation of a polaronic state and to charac- 
terize the static and dynamical properties of such a state. 
We will only be concerned with the Hubbard-Holstein 
and the t J-Holstein models, as minimal systems with lo- 
cal electron-electron (e-e) and electron-phonon (e-ph) in- 
teractions. In order to focus on the physics of unbound 
polarons (as opposed to bipolarons), we will always con- 
sider the limit of strong local e-e repulsion to prevent the 
formation of a local bipolaron |13| . The short-range char- 
acter of the bare interactions considered in this model 
greatly simplifies the numerical analysis on finite clusters. 



Furthermore, the strong coupling limit of the e-ph in- 
teraction gives rise to small (single-site) polarons, which 
are individually rather well understood, thus providing 
a simple limit of the model. Since we focus on metallic 
states, we purposely avoid specific conditions (like, e.g., 
quarter-filling) leading to ordered insulating states. Fi- 
nally we confine our analysis to the one-dimensional case. 
This choice is an inescapable consequence of the small- 
ness of our numerical clusters, but also allows for useful 
comparisons with the wealth of well-established physical 
results available in d = 1 . 
The model reads 



^ = -* E {^^,^^1^- + ^A + ^E^T^^.TcliC, 



^i,i 



{ij),<^ 



dYl [4<TCi,<T - (4^c,,^)] (fli + a^j +uJo^a\a,. (1) 



We use units such that the lattice spacing a — 1 and 
also Ti ^ c ^ \. To make our analysis more complete, 
we also investigate the iJ-Holstein model, where the 
U term in (|^) is replaced by an Heisenberg interaction 



K7 = JEfe)(s 



s,- 



iUj), with the additional con- 



straint of no double occupancy ^^ 



< 1. 



II. THE SMALL POLARON CRITERION 

A polaronic state can be characterized as a bound state 
of electrons and phonons, in which the electronic mo- 
tion is accompanied by a significant lattice displacement. 



The strong coupling between the electron and the lattice 
strongly suppresses the electronic mobility. The carriers 
are therefore self-trapped in the potential well that they 
generated. The polaronic state is therefore characterized 
by two conditions: (a) the energy gain associated with 
the self-trapping must be larger than the loss of kinetic 
energy, and (b) the local lattice displacement must be 
sizeable, in order to significantly reduce the electronic 
hopping amplitude. 

In previous analyses in the single-particle case of the 
Holstein model |g,|l2|,|[, emphasis was put on the differ- 
ence in the criteria for polaron formation in the adiabatic 
(luq < t ) or antiadiabatic (luq > t) regimes. In the adia- 
batic regime the crossover to polarons is dictated by the 
condition (a), since in this limit, as soon as a bound state 
is energetically favorable, the electron mobility is auto- 
matically reduced due to the large mass of the lattice. In 
the single particle case the energy of the strong-coupling 
polaronic state is Epoi — ^9^ /^o, while the free electron 
energy Eei = ^2dt. Then the condition for single polaron 
formation in the nearly adiabatic {too/t < 1) case reads 
(see, e.g. Ref. |). 



A = g^/{2dtu;o) > 1. 



(2) 



In the antiadiabatic regime the crossover for the single- 
particle case is instead given by the condition (b), which 
in the Holstein model is expressed as a = g/ujQ > 1. 

In the presence of many carriers and e-e correlations 
the electronic energies and the e-ph coupling are affected, 
so that the criterion for polaron formation in the adia- 
batic regime (which involves a condition on the energies) 
will be different in the many particle case with respect 
to the single particle one. On the other hand condition 
(b) ruling polaron formation in the antiadiabatic regime 
is not substantially changed by increasing the number of 
particles. In particular the slow dynamics of the elec- 
tronic degrees of freedom when luq ^ t strongly sup- 
presses electronic screening processes. Therefore in the 
antiadiabatic regime the single-particle condition a = 
g/ujo > 1 of large ion displacement is expected to hold 
also in the many particle case. 

For these reasons we mostly investigate here the phys- 
ically relevant case of the adiabatic regime where the 
phonon frequency ujo is smaller than the electronic en- 
ergy scale t (the typical electronic energy scale stays t 
even in the presence of correlations). The same analysis 
has also been applied to the antiadiabatic case, where 
it is confirmed that the single-particle condition a > 1 
holds in the many-particle case as well for all the con- 
sidered values of parameters and fillings. Therefore we 
will only present results for the adiabatic regime. In this 
case the energetic balance condition rules the crossover. 
A polaronic bound state can therefore be realized only 
if the energy of such a state Spoi is lower than the en- 
ergy of the electrons in the absence of electron-phonon 
interaction S^i ■ 



It is quite natural to generalize the condition (^ to 
the many-particle case by comparing the energy of the 
strongly coupled e-ph system, Spoi, where the polarons 
are strongly localized and the kinetic energy is negligi- 
ble, with the energy of the purely electronic system (i.e. 
of the simple Hubbard or tJ models) . In other words, 
we identify the value of the coupling for which the po- 
laronic crossover occurs, g = gc, with the c-ph coupling 
above which the energy of the strongly polaronic state is 
lower than the energy of the purely electronic state. In 
this framework it is useful to generalize the definition of 
A by defining the quantity 



A = £pol{g)/Sel 



(3) 



that reduces to A for a single particle, and show that the 
criterion for small polaron formation in the general case 
of finite densities and in the presence of e-e correlation 
is A > 1 . Of course this criterion can only be significant 
provided the crossover between the purely electronic and 
the polaronic regimes is sufhciently sharp, as is indeed 
the case in the adiabatic regime |g] . 

The criterion A = Spoi{g)/£ei > 1 would be of httle 
use if simple expressions in terms of the bare parame- 
ters or at least simple estimates of Spot (g) and £ei were 
not available. Therefore, we now turn to the explicit 
evaluation of Spoi (g) and £ei ■ Unfortunately, while some 
simple arguments allow an easy estimation of the energy 
in the strongly polaronic regime, the knowledge of the 
purely electronic energy (i.e. of the Hubbard model) is 
a much more difficult task. However, we point out that 
our goal is not to investigate the Hubbard model by itself, 
but rather to analyze the many-body effects of interac- 
tion on polaron formation, once a full knowledge of the 
purely electronic system (no matter how complicated) is 
assumed. 

In order to evaluate the strong coupling energy £poi{g) 
we start from the atomic limit of zero hopping i = 0, in 
which the Hubbard-Holstcin model can be exactly solved. 
In this case a Lang-Firsov canonical transformation, al- 
lows the elimination of the linear e-ph coupling by in- 
troducing a shift a = g/(-uo of the ionic equilibrium po- 
sition. An effective non-retarded attraction between the 
particles arises Uattr = —2g'^/uJo, which must be com- 
pensated by a larger repulsion U to avoid (local) bipo- 
laron formation. In this atomic limit each site is inde- 
pendent and each electron individually shifts the ionic 
site on which it resides forming a strong-coupling po- 
laron of energy Si^poi — —g'^/ajQ. An opposite shift, 
but with identical energy gain, is induced by holes, so 
that, for Ns sites at average filling n the total energy is 
^poi^ = —Nsn{l — n)g'^/u}o. As soon as a finite hopping is 
introduced, a coherent kinetic energy contribution arises, 
which, even in the absence of e-e interaction is exponen- 
tially small [t* = iexp(— a^)] and can safely be neglected 
in our estimate of the energy. However, more relevant 
additional corrections arise, which are due to incoher- 



ent virtual hopping processes of the electrons to neigh- 
boring sites. These processes naturally depend on the 
occupancy of the neighboring sites. If the site is unoc- 
cupied, the electron can hop onto it without paying any 
electronic repulsion U . However, this hopping occurs in- 
stantaneously without allowing the lattice time to relax. 
Therefore an intermediate virtual state of higher energy 
proportional to g'^/ujQ '^ t is reached, before a second 
hop restores the initial configuration. In the limit of very 
strong e-ph coupling and accounting for the probability 
for a singly occupied site to have at least one empty site 
nearby, the virtual double-hopping process lowers the to- 
tal energy by an amount of the order p4| 



fi/A « -2iV,n(l - n)f/{g^/LOo) 



(4) 



On the other hand a second incoherent correction to the 
polaronic energy occurs when the neighboring site is al- 
ready occupied by an electron with opposite spin (for 
parallel spins this process is forbidden by the Pauli prin- 
ciple). In this case a virtual state is reached where an 
e-e interaction is also present giving rise to a superex- 
change coupling. In the case of finite e-e repulsion U, 
this term provides an additional incoherent contribution 
to the energy of the polaronic state. In a recent work 
[ jl5[ we showed that the e-ph coupling dresses the effec- 
tive magnetic coupling J leading to an increase of its 
value. However, for a Holstein-like e-ph coupling no cor- 
rections to the purely electronic J arise up to order g^ /U. 
For the sake of simplicity we therefore estimate this con- 
tribution £'^i by using the value for ,9 = calculated by 
exact diagonalization. 

As a result, the energy of the polaronic state can be 
estimated as 



£. 



pol 



-N,n{l - n)gVwo + £i/x + £i 



(5) 



It is worth noting that in the above expression the bare 
e-ph coupling and phonon frequency appear. This is be- 
cause in the strong-coupling regime, where the coherent 
motion of the polarons is suppressed, no screening oc- 
curs due to the very massive carriers. This situation is 
different from the case of a Fermi liquid, where the pres- 
ence of the e-e interaction leads to a dressing of the e-ph 
coupling [|l6| and is complementary to the case of an in- 
teracting electron system near the Mott-Hubbard transi- 
tion (or near the Luther-Emery point for one-dimensional 
systems), where Umklapp processes render the e-ph in- 
teractions ineffective: If a system is localized by the e-e 
interaction, this latter screens out the effects of the e- 
ph coupling and, conversely, if a system is localized by 
strong e-ph interactions (as in the present case), the e-e 
interactions have no effect on the phonon parameters. 

The evaluation of the purely electronic energy E^i of 
the Hubbard model is not possible in general unless one 
resorts to numerical methods. However, for the sake 
of simplicity, we consider a one-dimensional system for 



which we can evaluate the energy in a simple way us- 
ing well-known results. In the U = oo limit an exact 
mapping exists between the Hubbard model and a sys- 
tem of free spinless fermions [[17f , and in the more general 
case of finite U the Bethe ansatz solution of the Hubbard 
model provides the exact results (see, e.g., Rcf. IQ). To 
make our calculation similar for the Hubbard and the tJ 
model, the alternative we chose is to directly carry out 
a numerical evaluation of the ground state energy of the 
purely electronic models on a one-dimensional cluster as 
a function of t/ or J and n. 

We are now in a position to check the validity of the 
criterion Ac = £poi{gc)/£ei '^ 1 by means of exact diag- 
onalization of finite clusters. Due to the infinite set of 
accessible phonon states on each site, we need to trun- 
cate the phonon Hilbert space by allowing for a finite 
maximum number NK""^ of phonons on each site. The 
crossover between quasi-free electrons and a polaronic 
state is signaled by rather abrupt changes in the behav- 
ior of most physical quantities; many quantities can be 
used to identify the crossover coupling, some examples of 
which are the electronic kinetic energy or effective mass, 
the electron-lattice correlation function or the average 
number of phonons per site, and so on. We identify the 
crossover value g — gc separating the two regimes with 
the value of g for which the slope of the average num- 
ber of phonons per site (riphig)) i^ maximum p(| . We 
always limited the calculation to values of g such that 
{''T'phig)) < ^ph°'^ and checked that the result was well 
converged by changing N^l"-^ . 

We performed calculation for various fillings in 5 and 
6 sites clusers for oji^/t = 0.2. In table I we report, be- 
sides the electronic energy and the magnetic energy, the 
crossover values of A and A obtained using the numeri- 
cally calculated values of gc- 

As far as A is concerned, we report results obtained 
both neglecting the incoherent 1/A contribution Eiix 
(sixth column, Ac) and considering such contribution (sev- 
enth column, A^) 

Note that for [/ = oo an additional symmetry ap- 
pears around quarter-filling in = f /2), making the filling 
n = f /5 equivalent to n = 4/5 and n — 2/5 equivalent to 
n = 3/5. 

It is evident from Table I that the criterion Ac « f is 
quite well satisfied, even neglecting the incoherent term, 
whereas the values of Ac are significantly different for 
different parameters. This result is particularly remark- 
able since the energies £poi(gc) and £ei entering in the 
numerator and the denominator of A respectively, vary 
significantly with filling and Coulomb repulsion U (or 
J). Various observations are in order. First of all it is 
apparent that the calculation of Ac without including the 
incoherent contribution £i/x to £poi gives quite homoge- 
neous values of Ac ~ f (within 10%) for systems with 
different values of U or J, but similar filling. Larger dis- 
crepancies (but always smaller than 20%) exist between 



different classes of fillings. We notice that the values 
of Ac for finite U or J are slightly underestimated be- 
cause we estimated the magnetic energy in Spoi taking 
the purely electronic contribution, which is known to be 
lower than in the presence of e-ph coupling |]l5| . This un- 
derestimation is proportionally more crude for larger J 
or smaller U values. Correcting for this effect should pro- 
duce an even more uniform value for Ac- The magnetic 
energy in the strongly polaronic regime is also different 
from £^i because the reduced itinerancy of the carriers 
obviously affects the spin correlations. Therefore a fur- 
ther refinement of the estimate of Spoi could be obtained 
by calculating the magnetic energy in the polaronic state 
taking E^^ from the purely electronic model with i = 0. 
Also in this case a larger magnetic coutribution would 
be found (when the carriers are localized the magnetic 
correlations are obviously stronger) and this would go in 
the direction of rendering the values for Ac more uniform. 
An additional important ingredient in the estimate of 
Ac turns out to be the contribution Sin. of the incoherent 
hopping processes to Spot [see Eq. (||)]. The estimated 
filling dependent contribution ^) leads to a sizable mod- 
ification of the criterion from Ac '^ 1 to Ac ~ 1.2. How- 
ever, despite this purely quantitative change, which was 
not considered in our previous single-polaron work [||, 
and which arises from a simple refinement of the esti- 
mate of the polaron energy, the physical meaning of the 
criterion stays the same. As a result of this latter filling- 
dependent improvement in the estimate of the polaron 
energy, the crossover values Ac's reported in the last col- 
umn of Table I turn out to be more homogeneous with 
varying filling within each class of interaction values and 
the criterion A'c « 1.2 seems to be generically satisfied 
within a few percent. 



III. DYNAMICAL PROPERTIES 



In the previous section, we characterized the crossover 
from weak to strong coupling regimes through the anal- 
ysis of static properties, namely the average number 
of phonons. In this section we consider the effects of 
this polaron-formation crossover on dynamical properties 
such as the spectral densities A{uj) and the dynamical 
conductivity cr(uj). Here we report on this investigation 
for the specific but typical case of the Hubbard-Holstein 
model with U — lOt. 
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A. Spectral Function 



In this section we study the spectral density associated 
with the injection of an electron (inverse photoemission. 



K')). 



(6) 



and the corresponding quantity for the emission of an 
electron (photoemission, PE) 



A-{uj) 
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In Fig.|] we show the PE (solid line) and IPE (dashed 
line) spectral densities for a 5-site cluster and four elec- 
trons with U = lot, for different values of g ranging from 
g = Oto g = 0.9t. 

A recent fc-resolved spectral analysis in the two- 
dimensional tJ model with two holes [gl|, shows a 
markedly different behavior between the incoherent (lo- 
cal) excitations at high energies, not substantially af- 
fected by e-ph coupling, and the dispersive "quasipar- 
ticle" states near the Fermi level, which become very 
massive with increasing g |22[ . On the contrary this 
behavior is not present in our one-dimensional cluster, 
where the fact that electrons are composites of holons 
and spinous leads to all electronic excitations in the up- 
per and lower Hubbard bands being incoherent. These 
are only marginally affected by the e-ph interaction and 
mainly display a broadening, which is particularly evi- 
dent in the upper Hubbard band. A closer inspection 
also allows one to detect the effects of multiphonon ex- 
citations dressing the purely electronic states by forming 
"shoulders" or "combs" with excitation energies spaced 
hy ujQ. 

The presence of a phonon-induced retarded local at- 
traction also reduces the local instantaneous Hubbard 
repulsion U. For the moderate g = 0.3t and g = 0.6t, 
this attractive effect is not large enough (< 2g'^/ujo < 2) 
to be clearly visible in a marked reduction of the gap 
between the lower and upper Hubbard bands. Rather 
a small leaking of spectral weight at the edges of the 
gap results from the much more evident broadening men- 
tioned above. The effective phonon-mediated attraction 
increases the amount of doubly occupied sites in the 
ground state thus modifying the relative weight of the 
spectral densities related to the injection (inverse pho- 
toemission, IPE) or to the emission (photoemission, PE) 
of electrons, represented by the dashed and solid lines 
in Fig. |l] respectively. For the large value oi U = lOt 
and for g = 0.3t and g = 0.6i, which are substantially 
smaller than those leading to bipolaron formation, the 
amount of doubly occupied sites in the ground state is 
always small and its small changes do not significantly af- 
fect the relative weight of the photoemission and inverse 
photoemission spectral features. As a consequence one 
observes an approximate, but quite well-satisfied conser- 
vation of the spectral weight in the PE and IPE parts 



of the spectrum separately for moderate values of g. As 
the e-ph couphng is further increased to g = 0.9t the en- 
ergy gap between the upper and the lower Hubbard band 
is substantially reduced and the weight conservation no 
longer holds separately for the PE and IPE parts of the 
spectrum. 

On the other hand, for the largest values of g we ob- 
serve a substantial transfer of spectral weight from the 
upper Hubbard band, which is broadened and lowers its 
overall intensity, to the lower Hubbard band. This lat- 
ter acquires weight essentially because the possibility of 
exciting phonons allows more states to be reached by in- 
jecting one electron in the system by conserving the mo- 
mentum and without increasing the number of doubly 
occupied sites in the final state. 



B. Optical Conductivity 

The real part of the conductivity for a one-dimensional 
tight-binding model at zero temperature may be ex- 
pressed in terms of the Kubo formula 



a{uj) = D6{uj) + 






K\3\cPo)rS{uj-{E„-E^')), (8) 



where J is the current operator. The finite frequency 
part (second term) can be straightforwardly calculated 
by means of the Lanczos algorithm. The coefficient of 
the zero-frequency delta function contribution D is usu- 
ally called the Drude weight and can be evaluated by 
combining the above with the well known sumrule for the 
total conductivity in terms of the kinetic energy (Ht), 



TTC 

a{uj)duj = ^{Ht). 



(9) 



In Fig. we show the behavior of the finite-frequency 
part of the dynamical conductivity for the same param- 
eters as in Fig. |l|. Again the most apparent effect of the 
coupling with the phonons is a substantial broadening of 
the high-energy excitations. A second remarkable feature 
is represented by the marked increase of spectral weight 
at low frequencies. As in the case of the single-particle 
spectral density, this occurs because the possibility of 
phonon excitations accompanying electron particle-hole 
excitations allows for a large number of intraband transi- 
tions (i.e. between states with the same number of doubly 
occupied sites). This effect is made particularly apparent 
by the comparison with the .g = case (panel (a) in Fig. 
^ , where the small number of sites and of holes strongly 
suppresses the weight at frequencies of order t |2^] . 

As far as the zero-frequency contribution is concerned, 
although quasiparticle transport is absent in d = 1, there 
is still coherent charge (holon) transport. At T=0 this 
leads to the presence of a (5-like Drude conductivity at 
zero frequency. It is noticeable that this ideal conduc- 
tance behavior persists at finite values of g at least at 



T = despite the non-integrable character of the model 
with 5 7^ 0. 

The suppression of coherent charge transport via 
phonon-dressing of the holon excitations leads to a strong 
decrease of the Drude spectral weight in the dynamical 
conductivity. As shown in Fig. y, the Drude weight re- 
lated to absorption from coherent excitations aX lo — Q 
decreases slowly in the weak-coupling regime (A < 1) 
whereas it is strongly suppressed for larger e-ph cou- 
pling when the polaronic crossover sets in. In the lat- 
ter case the rapid decrease is consistent with the usual 
exponential decrease of the charge-carrier band-width 

V -^ V exp —a^. The only difference is now that in rf = 1, 

V no longer is to be interpreted as the Fermi velocity of 
quasiparticles, but is related to the holon dispersion. On 
the other hand, by increasing g, spectral weight is trans- 
ferred to higher energies, as indicated by the increase of 
the finite frequency weight in Fig. 0. 

Although we cannot easily access the very strong cou- 
pling regime, we also checked in a few cases that for 
stronger couplings than those reported in the figure, this 
incoherent high-energy part of the absorption also de- 
creases, but more slowly, according to the expected 1/g 
behavior. The overall decrease of the total absorption is 
a natural consequence of the sum rule \^ relating the fre- 
quency integrated dynamical conductivity to the average 
kinetic energy. 



IV. CONCLUSIONS 

In this work, we have investigated the electron-polaron 
crossover in the many-particle case. We mainly consid- 
ered the adiabatic regime, in which the condition for po- 
laron formation is modified with respect to the single- 
particle case. The main result is that a clear criterion 
has been identified determining the strength of e-ph cou- 
pling leading to polaron formation. Remarkably, a simple 
physical interpretation of the criterion is still possible in 
the many-particle case provided that the single electron 
energy is replaced by the many-body purely electronic 
energy and the magnetic energy contribution is consid- 
ered in the polaronic phase. 

The effects of the intermediate polaronic regime have 
also been investigated on the dynamical properties. In 
particular, both in the single-particle spectral density 
and in the optical conductivity, we find that the electron- 
to-polaron crossover leads to a broadening of the high 
energy [uj ^ U) features. At the same time the in- 
crease of momentum conserving processes due to the 
mixed electron-multiphonon character of the excitations 
substantially increases the weight of intraband [uj ~ t) 
processes. 

The above scenario was investigated by numerical ex- 
act diagonalization of a 5-site cluster. This leads to 
several limitations. First of all the presence of finite 
size effects in such a small system forces the analysis 



to be limited to short-range models. Even the analy- 
sis of models with a nearest-neighbor intersite e-ph cou- 
pling like the Su-Schrieffer-Heeger model lead to sub- 
stantially larger finite-size effects. Another quite in- 
teresting problem, which we could not address within 
our small-cluster analysis is related to the competi- 
tion with other phonon-drivcn instabilities like charge- 
density- waves (CDW), phase separation or incommensu- 
rate stripe formation. Even in the case of these latter two 
instabilities, which are not Fermi-surface instabilities, the 
small size of our system does not allow for a systematic 
comparison with the polaron formation within our nu- 
meric analysis. On the other hand, some information 
can be gained by comparing our results with different 
analyses. In particular as far as phase separation and 
(with the inclusion of long-range Coulomb interactions) 
stripe formation, the analyses of Refs. P4 16 1 indicate 



that they occur for e-ph couplings that are smaller than 



those leading to polaron formation. The approximate 
character of these analyses, however, still leaves room for 
futher work in this direction. As far as commensurate 
CDW formation is concerned, the recent density-matrix 
renormalization group analysis of Ref. P5[ in the infinite- 
U Hubbard-Holstein model a,t n — 1/4 indicates that the 
CDW instability competes strongly with polaron forma- 
tion. However, in the adiabatic regime the CDW and 
the polaron instability occur nearby in the phase dia- 
gram (whereas polarons are favored in the antiadiabatic 
regime), so it is likely that the polaronic instability will 
win for fillings far from 1/4, which are less favorable for 
CDW. 

Part of this work has been carried out with the fi- 
nancial support of the INFM, PRA HTSC 1996. The 
Authors would like to acknowledge Prof. S. Ciuchi, C. 
Castellani, C. Di Castro and D. Feinberg for many fruit- 
ful discussions. 
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FIG. 1. Spectral density (solid line for PE, dashed line for IPE) for U = lOi and uo/t = 0.4, for various values of g and 
n = 4/5. 
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FIG. 3. Integrated optical conductivity (Total Sumrule), Drude weight and finite frequency weight as functions of g for 
V = lot, LJo/t = 0.4 and n = 4/5 
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TABLE I. Filling (first column) and interaction (second column) dependencies of the total purely electronic energy (third 
column), of the purely magnetic energy (fourth column), of Ac found with the maximum derivative of (riphig)) (fifth column), 
and of the Ac's calculated without{\c) and with (A'cthe £i/x correction in £poi for (jJo/t = 0.2 
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